Senescence defines a distinct subset of myofibroblasts that orchestrates immunosuppression in pancreatic cancer.

PDAC therapeutic resistance is largely attributed to a unique tumor microenvironment embedded with an abundance of cancer associated fibroblasts (CAFs). Distinct CAF populations were recently identified, but the phenotypic drivers and specific impact of CAF heterogeneity remain unclear. In this study, we identify a subpopulation of senescent myofibroblastic CAFs (SenCAFs) in mouse and human PDAC. These SenCAFs are a phenotypically distinct subset of myofibroblastic CAFs that localize near tumor ducts and accumulate with PDAC progression. To assess the impact of endogenous SenCAFs in PDAC, we employed a LSL-KRASG12D;p53flox;p48-CRE;INK-ATTAC (KPPC-IA) mouse model of spontaneous PDAC with inducible senescent cell depletion. Depletion of senescent stromal cells in genetic and pharmacologic PDAC models relieved immune suppression by macrophages, delayed tumor progression and increased responsiveness to chemotherapy. Collectively, our findings demonstrate that SenCAFs promote PDAC progression and immune cell dysfunction.

Senescence Defines a Distinct Subset of Myofibroblasts That Orchestrates Immunosuppression in Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) therapeutic resistance is largely attributed to a unique tumor microenvironment embedded with an abundance of cancer-associated fibroblasts (CAF). Distinct CAF populations were recently identified, but the phenotypic drivers and specific impact of CAF heterogeneity remain unclear. In this study, we identify a subpopulation of senescent myofibroblastic CAFs (SenCAF) in mouse and human PDAC. These SenCAFs are a phenotypically distinct subset of myofibroblastic CAFs that localize near tumor ducts and accumulate with PDAC progression. To assess the impact of endogenous SenCAFs in PDAC, we used an LSL-KRASG12D;p53flox;p48-CRE;INK-ATTAC (KPPC-IA) mouse model of spontaneous PDAC with inducible senescent cell depletion. Depletion of senescent stromal cells in genetic and pharmacologic PDAC models relieved immune suppression by macrophages, delayed tumor progression, and increased responsiveness to chemotherapy. Collectively, our findings demonstrate that SenCAFs promote PDAC progression and immune cell dysfunction. SIGNIFICANCE: CAF heterogeneity in PDAC remains poorly understood. In this study, we identify a novel subpopulation of senescent CAFs that promotes PDAC progression and immunosuppression. Targeting CAF senescence in combination therapies could increase tumor vulnerability to chemo- or immunotherapy. See related article by Ye et al.

Optimal cell traction forces in a generalized motor-clutch model.

Cells exert forces on mechanically compliant environments to sense stiffness, migrate, and remodel tissue. Cells can sense environmental stiffness via myosin-generated pulling forces acting on F-actin, which is in turn mechanically coupled to the environment via adhesive proteins, akin to a clutch in a drivetrain. In this "motor-clutch" framework, the force transmitted depends on the complex interplay of motor, clutch, and environmental properties. Previous mean-field analysis of the motor-clutch model identified the conditions for optimal stiffness for maximal force transmission via a dimensionless number that combines motor-clutch parameters. However, in this and other previous mean-field analyses, the motor-clutch system is assumed to have balanced motors and clutches and did not consider force-dependent clutch reinforcement and catch bond behavior. Here, we generalize the motor-clutch analytical framework to include imbalanced motor-clutch regimes, with clutch reinforcement and catch bonding, and investigate optimality with respect to all parameters. We found that traction force is strongly influenced by clutch stiffness, and we discovered an optimal clutch stiffness that maximizes traction force, suggesting that cells could tune their clutch mechanical properties to perform a specific function. The results provide guidance for maximizing the accuracy of cell-generated force measurements via molecular tension sensors by designing their mechanosensitive linker peptide to be as stiff as possible. In addition, we found that, on rigid substrates, the mean-field analysis identifies optimal motor properties, suggesting that cells could regulate their myosin repertoire and activity to maximize force transmission. Finally, we found that clutch reinforcement shifts the optimum substrate stiffness to larger values, whereas the optimum substrate stiffness is insensitive to clutch catch bond properties. Overall, our work reveals novel features of the motor-clutch model that can affect the design of molecular tension sensors and provide a generalized analytical framework for predicting and controlling cell adhesion and migration in immunotherapy and cancer.

Protocol to culture and image pancreatic ductal adenocarcinoma tumor slices to study T cell migration.

Here, we describe a protocol for culture and live cell imaging of tumor slices. This approach studies carcinoma and immune cell dynamics in complex tumor microenvironments (TME) with nonlinear optical imaging platforms. Using a tumor-bearing mouse model of pancreatic ductal adenocarcinoma (PDA), we detail steps to isolate, activate, and label CD8+ T lymphocytes and later introduce them to live murine PDA tumor slice explants. The techniques described in this protocol can improve our understanding of cell migration in complex microenvironments ex vivo. For complete details on the use and execution of this protocol, please refer to Tabdanov et al. (2021).1.

Differential contractility regulates cancer stem cell migration.

Cancer stem cells (CSCs) are known to have a high capacity for tumor initiation and the formation of metastases. We have previously shown that in collagen constructs mimetic of aligned extracellular matrix architectures observed in carcinomas, breast CSCs demonstrate enhanced directional and total motility compared with more differentiated carcinoma populations. Here, we show that CSCs maintain increased motility in diverse environments including on 2D elastic polyacrylamide gels of various stiffness, 3D randomly oriented collagen matrices, and ectopic cerebral slices representative of a common metastatic site. A consistent twofold increase of CSC motility across platforms suggests a general shift in cell migration mechanics between well-differentiated carcinoma cells and their stem-like counterparts. To further elucidate the source of differences in migration, we demonstrate that CSCs are less contractile than the whole population (WP) and develop fewer and smaller focal adhesions and show that enhanced CSC migration can be tuned via contractile forces. The WP can be shifted to a CSC-like migratory phenotype using partial myosin II inhibition. Inversely, CSCs can be shifted to a less migratory WP-like phenotype using microtubule-destabilizing drugs that increase contractility or by directly enhancing contractile forces. This work begins to reveal the mechanistic differences driving CSC migration and raises important implications regarding the potentially disparate effects of microtubule-targeting agents on the motility of different cell populations.

Targeting TNF-α-producing macrophages activates antitumor immunity in pancreatic cancer via IL-33 signaling.

Pancreatic ductal adenocarcinoma (PDA) remains resistant to immune therapies, largely owing to robustly fibrotic and immunosuppressive tumor microenvironments. It has been postulated that excessive accumulation of immunosuppressive myeloid cells influences immunotherapy resistance, and recent studies targeting macrophages in combination with checkpoint blockade have demonstrated promising preclinical results. Yet our understanding of tumor-associated macrophage (TAM) function, complexity, and diversity in PDA remains limited. Our analysis reveals significant macrophage heterogeneity, with bone marrow-derived monocytes serving as the primary source for immunosuppressive TAMs. These cells also serve as a primary source of TNF-α, which suppresses expression of the alarmin IL-33 in carcinoma cells. Deletion of Ccr2 in genetically engineered mice decreased monocyte recruitment, resulting in profoundly decreased TNF-α and increased IL-33 expression, decreased metastasis, and increased survival. Moreover, intervention studies targeting CCR2 with a new orthosteric inhibitor (CCX598) rendered PDA susceptible to checkpoint blockade, resulting in reduced metastatic burden and increased survival. Our data indicate that this shift in antitumor immunity is influenced by increased levels of IL-33, which increases dendritic cell and cytotoxic T cell activity. These data demonstrate that interventions to disrupt infiltration of immunosuppressive macrophages, or their signaling, have the potential to overcome barriers to effective immunotherapeutics for PDA.

Stromal architecture directs early dissemination in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDA) is an extremely metastatic and lethal disease. Here, in both murine and human PDA, we demonstrate that extracellular matrix architecture regulates cell extrusion and subsequent invasion from intact ductal structures through tumor-associated collagen signatures (TACS). This results in early dissemination from histologically premalignant lesions and continual invasion from well-differentiated disease, and it suggests TACS as a biomarker to aid in the pathologic assessment of early disease. Furthermore, we show that pancreatitis results in invasion-conducive architectures, thus priming the stroma prior to malignant disease. Analysis in potentially novel microfluidic-derived microtissues and in vivo demonstrates decreased extrusion and invasion following focal adhesion kinase (FAK) inhibition, consistent with decreased metastasis. Thus, data suggest that targeting FAK or strategies to reengineer and normalize tumor microenvironments may have roles not only in very early disease, but also for limiting continued dissemination from unresectable disease. Likewise, it may be beneficial to employ stroma-targeting strategies to resolve precursor diseases such as pancreatitis in order to remove stromal architectures that increase risk for early dissemination.

Aligned forces: Origins and mechanisms of cancer dissemination guided by extracellular matrix architecture.

Organized extracellular matrix (ECM), in the form of aligned architectures, is a critical mediator of directed cancer cell migration by contact guidance, leading to metastasis in solid tumors. Current models suggest anisotropic force generation through the engagement of key adhesion and cytoskeletal complexes drives contact-guided migration. Likewise, disrupting the balance between cell-cell and cell-ECM forces, driven by ECM engagement for cells at the tumor-stromal interface, initiates and drives local invasion. Furthermore, processes such as traction forces exerted by cancer and stromal cells, spontaneous reorientation of matrix-producing fibroblasts, and direct binding of ECM modifying proteins lead to the emergence of collagen alignment in tumors. Thus, as we obtain a deeper understanding of the origins of ECM alignment and the mechanisms by which it is maintained to direct invasion, we are poised to use the new paradigm of stroma-targeted therapies to disrupt this vital axis of disease progression in solid tumors.

Engineering T cells to enhance 3D migration through structurally and mechanically complex tumor microenvironments.

Defining the principles of T cell migration in structurally and mechanically complex tumor microenvironments is critical to understanding escape from antitumor immunity and optimizing T cell-related therapeutic strategies. Here, we engineered nanotextured elastic platforms to study and enhance T cell migration through complex microenvironments and define how the balance between contractility localization-dependent T cell phenotypes influences migration in response to tumor-mimetic structural and mechanical cues. Using these platforms, we characterize a mechanical optimum for migration that can be perturbed by manipulating an axis between microtubule stability and force generation. In 3D environments and live tumors, we demonstrate that microtubule instability, leading to increased Rho pathway-dependent cortical contractility, promotes migration whereas clinically used microtubule-stabilizing chemotherapies profoundly decrease effective migration. We show that rational manipulation of the microtubule-contractility axis, either pharmacologically or through genome engineering, results in engineered T cells that more effectively move through and interrogate 3D matrix and tumor volumes. Thus, engineering cells to better navigate through 3D microenvironments could be part of an effective strategy to enhance efficacy of immune therapeutics.

Engineering Elastic Nano- and Micro-Patterns and Textures for Directed Cell Motility.

We present a reproducible protocol for fabrication of polyacrylamide (PAA) hydrogel-based nano-patterns and nano-textures with a wide range of elastic rigidities to study fundamental cell behaviors, such as mechanosensitivity and motility. We explore the benefits of this protocol by successfully testing the compatibility of the PAA platforms with super-resolution microscopy, which is largely unavailable with platforms of nano-scale textures made from different polymers. We also utilized soft and rigid nano-textures to study the mechanosensing basis of T cell behavior and phenotype. For complete information on the generation and use of this protocol, please refer to Tabdanov et al. (2018b).

Loss of HIF1A From Pancreatic Cancer Cells Increases Expression of PPP1R1B and Degradation of p53 to Promote Invasion and Metastasis.

BACKGROUND & AIMS: Pancreatic ductal adenocarcinomas (PDACs) are hypovascular, resulting in the up-regulation of hypoxia inducible factor 1 alpha (HIF1A), which promotes the survival of cells under low-oxygen conditions. We studied the roles of HIF1A in the development of pancreatic tumors in mice. METHODS: We performed studies with KrasLSL-G12D/+;Trp53LSL-R172H/+;Pdx1-Cre (KPC) mice, KPC mice with labeled pancreatic epithelial cells (EKPC), and EKPC mice with pancreas-specific depletion of HIF1A. Pancreatic and other tissues were collected and analyzed by histology and immunohistochemistry. Cancer cells were cultured from PDACs from mice and analyzed in cell migration and invasion assays and by immunoblots, real-time polymerase chain reaction, and liquid chromatography-mass spectrometry. We performed studies with the human pancreatic cancer cell lines PATU-8988T, BxPC-3, PANC-1, and MiaPACA-2, which have no or low metastatic activity, and PATU-8988S, AsPC-1, SUIT-2 and Capan-1, which have high metastatic activity. Expression of genes was knocked down in primary cancer cells and pancreatic cancer cell lines by using small hairpin RNAs; cells were injected intravenously into immune-competent and NOD/SCID mice, and lung metastases were quantified. We compared levels of messenger RNAs in pancreatic tumors and normal pancreas in The Cancer Genome Atlas. RESULTS: EKPC mice with pancreas-specific deletion of HIF1A developed more advanced pancreatic neoplasias and PDACs with more invasion and metastasis, and had significantly shorter survival times, than EKPC mice. Pancreatic cancer cells from these tumors had higher invasive and metastatic activity in culture than cells from tumors of EKPC mice. HIF1A-knockout pancreatic cancer cells had increased expression of protein phosphatase 1 regulatory inhibitor subunit 1B (PPP1R1B). There was an inverse correlation between levels of HIF1A and PPP1R1B in human PDAC tumors; higher expression of PPP1R1B correlated with shorter survival times of patients. Metastatic human pancreatic cancer cell lines had increased levels of PPP1R1B and lower levels of HIF1A compared with nonmetastatic cancer cell lines; knockdown of PPP1R1B significantly reduced the ability of pancreatic cancer cells to form lung metastases in mice. PPP1R1B promoted degradation of p53 by stabilizing phosphorylation of MDM2 at Ser166. CONCLUSIONS: HIF1A can act a tumor suppressor by preventing the expression of PPP1R1B and subsequent degradation of the p53 protein in pancreatic cancer cells. Loss of HIF1A from pancreatic cancer cells increases their invasive and metastatic activity.

Fibrillar Collagen Quantification With Curvelet Transform Based Computational Methods.

Quantification of fibrillar collagen organization has given new insight into the possible role of collagen topology in many diseases and has also identified candidate image-based bio-markers in breast cancer and pancreatic cancer. We have been developing collagen quantification tools based on the curvelet transform (CT) algorithm and have demonstrated this to be a powerful multiscale image representation method due to its unique features in collagen image denoising and fiber edge enhancement. In this paper, we present our CT-based collagen quantification software platform with a focus on new features and also giving a detailed description of curvelet-based fiber representation. These new features include C++-based code optimization for fast individual fiber tracking, Java-based synthetic fiber generator module for method validation, automatic tumor boundary generation for fiber relative quantification, parallel computing for large-scale batch mode processing, region-of-interest analysis for user-specified quantification, and pre- and post-processing modules for individual fiber visualization. We present a validation of the tracking of individual fibers and fiber orientations by using synthesized fibers generated by the synthetic fiber generator. In addition, we provide a comparison of the fiber orientation calculation on pancreatic tissue images between our tool and three other quantitative approaches. Lastly, we demonstrate the use of our software tool for the automatic tumor boundary creation and the relative alignment quantification of collagen fibers in human breast cancer pathology images, as well as the alignment quantification of in vivo mouse xenograft breast cancer images.

Modeling distributed forces within cell adhesions of varying size on continuous substrates.

Cell migration and traction are essential to many biological phenomena, and one of their key features is sensitivity to substrate stiffness, which biophysical models, such as the motor-clutch model and the cell migration simulator can predict and explain. However, these models have not accounted for the finite size of adhesions, the spatial distribution of forces within adhesions. Here, we derive an expression that relates varying adhesion radius ( R) and spatial distribution of force within an adhesion (described by s) to the effective substrate stiffness ( κsub ), as a function of the Young's modulus of the substrate ( E Y ), which yields the relation, κsub=RsEY , for two-dimensional cell cultures. Experimentally, we found that a cone-shaped force distribution ( s = 1.05) can describe the observed displacements of hydrogels deformed by adherent U251 glioma cells. Also, we found that the experimentally observed adhesion radius increases linearly with the cell protrusion force, consistent with the predictions of the motor-clutch model with spatially distributed clutches. We also found that, theoretically, the influence of one protrusion on another through a continuous elastic environment is negligible. Overall, we conclude cells can potentially control their own interpretation of the mechanics of the environment by controlling adhesion size and spatial distribution of forces within an adhesion.

The role of nonmuscle myosin 2A and 2B in the regulation of mesenchymal cell contact guidance.

Contact guidance refers to the ability of cells to sense the geometrical features of the microenvironment and respond by changing their shape and adopting the appropriate orientation. Inhibition and ablation of nonmuscle myosin 2 (NM2) paralogues have demonstrated their importance for contact guidance. However, the specific roles of the NM2 paralogues have not been systematically studied. In this work we use micropatterned substrates to examine the roles of NM2A and NM2B and to elucidate the relationship of the microenvironment, actomyosin, and microtubules in contact guidance. We show that contact guidance is preserved following loss of NM2B and that expression of NM2A alone is sufficient to establish an appropriate orientation of the cells. Loss of NM2B and overexpression of NM2A result in a prominent cell polarization that is found to be linked to the increased alignment of microtubules with the actomyosin scaffold. Suppression of actomyosin with blebbistatin reduces cell polarity on a flat surface, but not on a surface with contact guidance cues. This indicates that the lost microtubule-actomyosin interactions are compensated for by microtubule-microenvironment interactions, which are sufficient to establish cell polarity through contact guidance.

Non-Invasive Monitoring of Stromal Biophysics with Targeted Depletion of Hyaluronan in Pancreatic Ductal Adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDA) is characterized by a pronounced fibroinflammatory stromal reaction consisting of inordinate levels of hyaluronan (HA), collagen, immune cells, and activated fibroblasts that work in concert to generate a robust physical barrier to the perfusion and diffusion of small molecule therapeutics. The targeted depletion of hyaluronan with a PEGylated recombinant human hyaluronidase (PEGPH20) lowers interstitial gel-fluid pressures and re-expands collapsed intratumoral vasculature, improving the delivery of concurrently administered agents. Here we report a non-invasive means of assessing biophysical responses to stromal intervention with quantitative multiparametric magnetic resonance imaging (MRI) at 14 Tesla (T). We found that spin-spin relaxation time T2 values and glycosaminoglycan chemical exchange saturation transfer (GagCEST) values decreased at 24 h, reflecting depletion of intratumoral HA content, and that these parameters recovered at 7 days concurrent with replenishment of intratumoral HA. This was also reflected in an increase in low-b apparent diffusion coefficient (ADC) at 24 h, consistent with improved tumor perfusion that again normalized at 7 days after treatment. Phantom imaging suggests that the GagCEST signal is driven by changes in HA versus other glycosaminoglycans. Thus, multiparametric magnetic resonance imaging (MRI) can be used as a non-invasive tool to assess therapeutic responses to targeted stromal therapy in PDA and likely other stroma-rich solid tumors that have high levels of hyaluronan and collagen.

Antifibrotic Therapy Disrupts Stromal Barriers and Modulates the Immune Landscape in Pancreatic Ductal Adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDA) remains one of the deadliest forms of cancer, in part, because it is largely refractory to current therapies. The failure of most standard therapies in PDA, as well as promising immune therapies, may be largely ascribed to highly unique and protective stromal microenvironments that present significant biophysical barriers to effective drug delivery, that are immunosuppressive, and that can limit the distribution and function of antitumor immune cells. Here, we utilized stromal reengineering to disrupt these barriers and move the stroma toward normalization using a potent antifibrotic agent, halofuginone. In an autochthonous genetically engineered mouse model of PDA, halofuginone disrupted physical barriers to effective drug distribution by decreasing fibroblast activation and reducing key extracellular matrix elements that drive stromal resistance. Concomitantly, halofuginone treatment altered the immune landscape in PDA, with greater immune infiltrate into regions of low hylauronan, which resulted in increased number and distribution of both classically activated inflammatory macrophages and cytotoxic T cells. In concert with a direct effect on carcinoma cells, this led to widespread intratumoral necrosis and reduced tumor volume. These data point to the multifunctional and critical role of the stroma in tumor protection and survival and demonstrate how compromising tumor integrity to move toward a more normal physiologic state through stroma-targeting therapy will likely be an instrumental component in treating PDA. SIGNIFICANCE: This work demonstrates how focused stromal re-engineering approaches to move toward normalization of the stroma disrupt physical barriers to effective drug delivery and promote antitumor immunity.See related commentary by Huang and Brekken, p. 328.

Bimodal sensing of guidance cues in mechanically distinct microenvironments.

Contact guidance due to extracellular matrix architecture is a key regulator of carcinoma invasion and metastasis, yet our understanding of how cells sense guidance cues is limited. Here, using a platform with variable stiffness that facilitates uniaxial or biaxial matrix cues, or competing E-cadherin adhesions, we demonstrate distinct mechanoresponsive behavior. Through disruption of traction forces, we observe a profound phenotypic shift towards a mode of dendritic protrusion and identify bimodal processes that govern guidance sensing. In contractile cells, guidance sensing is strongly dependent on formins and FAK signaling and can be perturbed by disrupting microtubule dynamics, while low traction conditions initiate fluidic-like dendritic protrusions that are dependent on Arp2/3. Concomitant disruption of these bimodal mechanisms completely abrogates the contact guidance response. Thus, guidance sensing in carcinoma cells depends on both environment architecture and mechanical properties and targeting the bimodal responses may provide a rational strategy for disrupting metastatic behavior.

Microtubule-Actomyosin Mechanical Cooperation during Contact Guidance Sensing.

Cancer cell migration through and away from tumors is driven in part by migration along aligned extracellular matrix, a process known as contact guidance (CG). To concurrently study the influence of architectural and mechanical regulators of CG sensing, we developed a set of CG platforms. Using flat and nanotextured substrates with variable architectures and stiffness, we show that CG sensing is regulated by substrate stiffness and define a mechanical role for microtubules and actomyosin-microtubule interactions during CG sensing. Furthermore, we show that Arp2/3-dependent lamellipodia dynamics can compete with aligned protrusions to diminish the CG response and define Arp2/3- and Formins-dependent actin architectures that regulate microtubule-dependent protrusions, which promote the CG response. Thus, our work represents a comprehensive examination of the physical mechanisms influencing CG sensing.

Cancer Stem Cell Migration in Three-Dimensional Aligned Collagen Matrices.

Cell migration is strongly influenced by the organization of the surrounding 3-D extracellular matrix. In particular, within fibrous solid tumors, carcinoma cell invasion may be directed by patterns of aligned collagen in the extra-epithelial space. Thus, studying the interactions of heterogeneous populations of cancer cells that include the stem/progenitor-like cancer stem cell subpopulation and aligned collagen networks is critical to our understanding of carcinoma dissemination. Here, we describe a robust method to generate aligned collagen matrices in vitro that mimic in vivo fiber organization. Subsequently, a protocol is presented for seeding aligned matrices with distinct carcinoma cell subpopulations and performing live cell imaging and quantitative analysis of cell migration. Together, the engineered constructs and the imaging techniques laid out here provide a platform to study cancer stem cell migration in 3-D anisotropic collagen with real-time visualization of cellular interactions with the fibrous matrix. © 2018 by John Wiley & Sons, Inc.

Dynamics of 3D carcinoma cell invasion into aligned collagen.

Carcinoma cells frequently expand and invade from a confined lesion, or multicellular clusters, into and through the stroma on the path to metastasis, often with an efficiency dictated by the architecture and composition of the microenvironment. Specifically, in desmoplastic carcinomas such as those of the breast, aligned collagen tracks provide contact guidance cues for directed cancer cell invasion. Yet, the evolving dynamics of this process of invasion remains poorly understood, in part due to difficulties in continuously capturing both spatial and temporal heterogeneity and progression to invasion in experimental systems. Therefore, to study the local invasion process from cell dense clusters into aligned collagen architectures found in solid tumors, we developed a novel engineered 3D invasion platform that integrates an aligned collagen matrix with a cell dense tumor-like plug. Using multiphoton microscopy and quantitative analysis of cell motility, we track the invasion of cancer cells from cell-dense bulk clusters into the pre-aligned 3D matrix, and define the temporal evolution of the advancing invasion fronts over several days. This enables us to identify and probe cell dynamics in key regions of interest: behind, at, and beyond the edge of the invading lesion at distinct time points. Analysis of single cell migration identifies significant spatial heterogeneity in migration behavior between cells in the highly cell-dense region behind the leading edge of the invasion front and cells at and beyond the leading edge. Moreover, temporal variations in motility and directionality are also observed between cells within the cell-dense tumor-like plug and the leading invasive edge as its boundary extends into the anisotropic collagen over time. Furthermore, experimental results combined with mathematical modeling demonstrate that in addition to contact guidance, physical crowding of cells is a key regulating factor orchestrating variability in single cell migration during invasion into anisotropic ECM. Thus, our novel platform enables us to capture spatio-temporal dynamics of cell behavior behind, at, and beyond the invasive front and reveals heterogeneous, local interactions that lead to the emergence and maintenance of the advancing front.